/// 公式为： 1 - (365 * 364 * ... (365 - N + 1) / 365 ^ N)
/// 上面公式在程序中的计算为： 1 - (365 / 365) * (364 / 365) * .. ((365 - N + 1) / 365)
pub fn new_birthday_probability(n: u32) -> f64 {
    assert!(n >= 2);
    let days_on_year = 365f64; // 这里没有考虑闰年
    let mut temp = days_on_year;
    let mut result = 1.;
    for _ in 0..n {
        result *= temp / days_on_year;
        temp -= 1.;
    }
    1. - result
}
